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Abstract 

MCSCF calculations are performed in order to determine the exchange coupling between the 
2p electrons of the N atom and the LUMOs of the fullerene cage in the case of mono- and tri- 
anions of NQCeo- The exchange coupling resulted by our calculations is large compared to the 
hyperfine interaction. The strong coupling can explain the missing EPR signal of the nitrogen in 
paramagnetic anions. 



PACS numbers: 31.15.Ar 



I. INTRODUCTION 



Since the discovery of the first endohedral fullerene a lot of interests have been attracted 
by this area of the nanotechnology. Many metal atom can be encapsulated by using discharge 
techniques or ion implantation. In all the cases the metal atom interacts strongly with the 
fullerene and acts as an electron donor occupying an 'off-centered' position inside the cage. 
In contrast the nitrogen in NtQiCeo is situated at the center of the molecule and retains its 
S=3/2 spin quartet atomic state^. This amazing property of the encapsulated N atom trig- 
gered several research on its possible application in quantum computing and spin labeling. 
Several publicationa^i^i^i^ studied the promise and limitations of using endohedral fuUerenes 
as quantum information carriers. Mehring et al.— recently pointed out experimentally the 
entanglement of the nuclear spin and the electronic spin of the encaged N atom. 

The changes of the characteristic EPR signal of the quartet electronic spin of the N atom 
makes it an ideal probe for monitoring chemical reactions of C60^. During the last decade 
a great deal of excitement has been brought by the discovery of the superconductivity of 
the alkali-doped fuUerenes. In this type of fullerene compounds the valence electrons of the 
ionized alkali atoms partially occupy the bands formed by the LUMOs of the Cgo molecules. 
The applicability of the quartet atomic state of the N atom as a spin label depends on the 
strength of the interaction between the 2p electrons of the N atom and the valence electron 
of the fullerene cage. An interaction which is small compared to the hyper-fine interaction, 
results in a line width effect of the EPR signal and the N^lCeo is a good candidate for a spin 
labeling agent. In the case of strong coupling the EPR signal of the system is completely 
changed and the lines corresponding to the valence electrons of the N atom are hard to 
identify in the signal of the paramagnetic system. 

The interaction between the 2p electrons of the N atom and the valence electrons of the 
Cgo can be described by a Heisenberg like effective Hamiltonian Hint = J^n^Cgo where 
Sat and Scgo denote the spin operator for the valence electrons of the N atom and the 
the Ceo, respectively, and J is the exchange coupling characterizing the strength of the 
interaction. The aim of the present paper is to determine theoretically the exchange coupling 
of the effective Hamiltonian. The exchange coupling has importance not only for EPR 
measurements but it plays essential role in the description of the transport through magnetic 
molecules^ which is particularly interesting from the point of view of spintronics. 
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II. COMPUTATIONAL DETAILS 



The calculations have been performed using the Gamess quantum chemical program 
package^. The proper description of the open-shell N@Cgo^ and N@CgQ^ anions requires multi- 
determinant wave functions. The restricted open-shell (ROHF) calculations in the Gamess 
package are accessible via the generalized valence bond (GVB) or the multi-configurational 
self-consistent field (MCSCF) methods using an appropriate active space. The energy of 
the anions of NQCeo with different multiplicity has been determined by means of CAS SCF 
calculations where the active space is confined to the 2p orbitals of the nitrogen atom and the 
three fold degenerate LUMOs of the fullerene molecule. For the clear interpretation of the 
results the excitations from the orbitals of the nitrogen to the LUMOs of the cage, and vice 
versa, were excluded from the active space applying the occupation restricted multiple active 
spacei^ technique. The MCSCF treatment of the open-shell systems using such a small active 
space is practically equivalent to the ROHF level of calculations. The calculation has been 
performed using split valence 631g basis on the carbon atoms. For the better description of 
the week interaction between the encapsulated atom and the fullerene molecule^i the basis 
on the N atom is extended by additional diffuse p orbitals and two d polarization functions 
(631-|-g(dd)). It is well known that in order to describe the electronic structure of negatively 
charged species application of diffuse basis functions is necessary. In our case the excess 
charge is distributed uniformly among the 60 carbon atoms and the lack of the diffuse basis 
on the carbon atoms does not affect dramatically our results. However, in order to check 
the sensitivity of the exchange coupling to the applied basis the calculations have also been 
performed with the Dunning's double zeta^^ and the split valence 631+g basis sets on the 
carbon atoms. The geometry of the N@CgQ^ molecule in the S=l state and the N@Cgo^ in 
the high spin S=3 state have been optimized at ROHF level and it is retained during the 
calculations of the energy of the systems with different multiplicity. 



III. RESULTS AND DISCUSSIONS 



It has been shown experimentally that in the highly reduced states of the N^Ceo the 
excess electrons occupy the LUMOs of the fullerene and the N atom inside the cage remains in 
spin quartet stated. In order to check the consistency of our calculations to the experimental 
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findings we performed a set of ROHF calculations on the mono- and tri-anions populating at 
first the 2p orbitals of the nitrogen and then populating the LUMOs of the Ceo- The results 
are summarized in Table [B The valence electrons of the nitrogen referred as N2p in Table [1] 
occupy the 7tiu orbitals of the endohedral complex between the Gh^ HOMO and 8tiu LUMO 
of the Geo in agreement with the result of rei^. Rather different value for the one-electron 
energy of the N2p orbitals is reported by Greer^^. This discrepancy is originated from the 
different treatment of the open-shell problem as it is discussed in ref.— . In the case of the 
mono-anion the energy of the two triplet states were compared while in the case of the triply 
ionized molecule the energy of the singlet state with fully occupied valence orbitals of N was 
compared to the high spin state of the NOGg^f . For both ions the system with intact N atom 
were energetically more favorable in agreement with the EPR measurements^. 

The interaction between the electrons of the nitrogen atom and the valence electrons on 
the Geo anion is described by a Heisenberg-like effective Hamiltonian: 

Hint = J^N^dio (1) 

where J is the coupling constant, Sn and Scgg are the spin of the nitrogen atom and the Geo 
anion, respectively. The square of the total spin operator = (S^v + Scgo)^ commutes with 
the Hamiltonian of the full system H = Hj^ + Hcgo + Hint and , consequently, its eigenvalue is 
a good quantum number. Expressing the interaction in terms of the spin of the subsystems 
and the spin of the whole molecule: 

Hint = - J {S^ -S% - S^gJ (2) 

the energy can be simply given as: 

Es = Eo + ^JSiS + l) , (3) 

where Eq denotes the energy of the separated systems and the subscript S indicates the 
explicit dependence of the energy on the multiplicity. Since our interaction Hamiltonian can 
describe only such processes in which Sn and Scgp are unchanged the excitations altering 
the spin of the subsystems, namely the hole and the particle are on different species, has to 
be excluded from the configuration space. In the case of N@GgQ^ Sjy = 3/2 and Scao = 1/2 
spanning an 8 dimensional direct product space. The total spin can have the values of 5 = 1 
or S* = 2 with the corresponding energies 

Es=i =Eo + J, Es=2 = Eo + 3J. (4) 
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Comparing the energy of the triplet and quintet state one can easily extract the exchange 
coupling as: 

J=^{Es=2-Es=i) (5) 

The results of the MCSCF calculations using 631g and DH basis are summarized in Table 
[m Although the application of the double zeta basis resulted in considerably deeper total 
energy the deviation of the exchange couplings is small. 

In the case of the triply ionized NtQCgo the valence electrons form a Scao =3/2 state on 
the LUMOs of the fuUerene molecule according to the Hund's rule. From the two quartet 
states, Sat, S^gg, four eigenstate of the operator can be constructed with the spin of 
S = 0,1,2,3, respectively. The corresponding energies as a function of S must be on a 
parabola according to Eq. [3l The results provided by the MCSCF calculations using three 
different basis sets are shown by Fig. [H The energies can nicely be fitted by the parabola 
given by Eq. [3l The exchange couplings obtained by using the split valence basis with 
and without diffuse p orbitals are practically the same. The inclusion of the diffuse basis 
functions on the carbon atoms resulted in negligible change. Although the magnitude of the 
exchange coupling corresponding to the double zeta basis is somewhat smaller than those 
provided by the split valence basis the agreement between them is satisfactory. 

Ferromagnetic exchange couplings between the 2p orbitals of the N atom and the valence 
electrons of the fuUerene molecule have been found in both anions. The exchange coupling 
of approximately 1 meV provided by our calculations for both systems is within the range of 
those found in organic ferromagnet&l^. This relatively strong coupling between the valence 
electrons of the nitrogen and the valence electrons of the fuUerene cage could be responsible 
for the disappearance of the nitrogen lines in the EPR spectrum of NQCgo anions with 
partially filled LUMOai^. 

IV. CONCLUSIONS 

In conclusion, ROHF and MCSF calculations have been performed on singly and triply 
ionized anions of NtQCgo in order to determine the effective exchange coupling between the 
valence electrons of the encapsulated N atom and the fullerene cage. In agreement with 
experiments we found that the excess electrons occupy the LUMOs of the fullerene molecule 
and the entrapped atom keeps its atomic character. The interaction between the valence 



5 



electrons of the N atom and the LUMOs of the Ceo can be well described by a Heisenberg 
like Hamiltonian. The size of the exchange couplings obtained by our calculations are much 
larger then the hyperfine interaction and can explain the results of EPR measurements on 
radical anions of NCijCeo- 

V. ACKNOWLEGMENTS 

This work is supported by the Hungarian National Science Foundation (contracts OTKA 
T038191 and T037856). 



1 T.Almeida Murphy, T. Pawlik, A. Weidinger, M. Hohne, R. Alcala, J.-M. Spaeth, Phys.Rev. 
Letter. 77 (1996) 1075 B. Pietzak, M. Waiblinger, T.Almeida Murphy, A. Weidinger, M. Hohne, 
E. Dietel, A. Hirsch, Chem. Phys. Letters 279 (1997) 259 

2 W. Harneit, Phys. Rev. A 65 (2002) 032322 

3 D. Suter and K. Lim, Phys. Rev. A 65 (2002) 052309 
J. Twamley Phys. Rev. A 67 (2003) 052318 

5 M. Feng and J. Twamley Phys. Rev. A 70 (2004) 032318 

6 M. Mehring, W. Scherer, and A. Weidinger, Phys. Rev. Lett. 93 (2004) 206603 

^ E. Dietel, A. Hirsch, B. Pietzak, M. Waiblinger, K. Lips, A. Weidinger, A. Gruss, K.P. Dinse, 

J. Am. Chem. Soc. 121 (1999) 2432 
^ F. Elste and C. Timm, Phys. Rev. B 71 (2005) 155403 

^ M.W.Schmidt, K.K.Baldridge, J.A.Boatz, S.T.Elbert, M.S.Gordon, J.H.Jensen, S.Koseki, 

N.Matsunaga, K.A.Nguyen, S.J.Su, T.L.Windus, M.Dupuis, 

J.A.Montgomery J.Comput.Chem. 14, 1347-1363(1993) 
° J.Ivanic J.Chem.Phys. 119 ((2003) 9364, 9377 

^ J.M. Park, P. Tarakeshwar, and K.S. Kim J. Chem. Phys.116 (2002) 10684 

^ T.H. Dunning, Jr., P.J. Hay Chapter 1 in "Methods of Electronic Structure Theory", H.F.Shaefer 

HI, Ed. Plenum Press, N.Y. 1977, pp 1-27. 
^ KP. Dinse, B. Godde, P. Jakes, M. Waiblinger, A. Weidinger, A. Hirsch, Abstr. Pap. - Am. 

Chem. Soc. (2001) 221st IEC-199. CODEN: ACSRAL ISSN: 0065-7727., P. Jakes, B. Godde, 



6 



M. Waiblinger, N. Weiden, K.P. Dinse, A. Weidinger, AIP Conference Proceedings 544 (2000) 
174 

D. A. Shultz,K.E. Vostrikova, S.H. Bodnar, Hyun-Joo Koo, Myung-Hwan Whangbo, M.L. Kirk, 

E. G. Depperman, and J.W. Kampf J. Am. Clicm. Soc. 125 (2003) 1607 
J. Lu, X. Zhang, X. Zhao, Chem. Phys. Lett. 312 (1999) 85 

J.C. Greer, Chem. Phys. Lett. 326 (2000) 567 

B.N. Plakhutin, N.N. Breslavskaya, E.V. Gorelik, A.V. Arbuznikov, Journal of Molecular Struc- 
ture: THEOCHEM 727 (2005) 149 



7 



Configuration E^ta; (Hartree) A E (eV) 



m 




N2/C6o8t?„ 
N2p3C6o8t{„ 


S = 
s = 


1 
1 


-2325.30365 
-2325.37155 


-1.84 


(a) 
(b) 


m 




N2/C6o8t?„ 


s = 
s = 



3 


-2324.74514 
-2325.10707 


-9.84 


(a) 
(b) 



TABLE I: Energies of N@Cqq and NQCgff with excess electron(s) occupying the 2p orbitals of 
the N atom (a) and the LUMOs of the Ceo niolecule (6). 



basis E5=i (Hartree) E5=2 (Hartree) J (meV) 
631g -2325.371558 -2325.371673 -1.56 
DH -2325.515025 -2325.515134 -1.49 



TABLE II: Energy of the N@Cqq resulted by MCSCF calculations using split valence (631g) and 
double zeta (DH) basis on the carbon atoms and the exchange coupling extracted from the energies. 



> 

E 



^ 



-1 ■ 
-2 ■ 
-3 ■ 
-4 ■ 
-5 ■ 
-6 ■ 
-7 ■ 
-8 ■ 



J631+g = -1-30 meV 
J63ig =-■'■31 "leV 



''DH 



-1.15 meV 



DH + 
631 g X 
631 +g M 



FIG. 1: Energies of N@CgQ^ corresponding to different multiplicity and the parabola fitted to the 



data points. The energy Eq independent of spin is subtracted. 
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